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ABSTRACT 



We present an improved version of SOAP (Boisse et al. 2012) named "SOAP-T" , which can generate the radial velocity variations 
and light-curves for systems consisting of a rotating spotted star with a transiting planet. This tool can be used to study the 
anomalies inside transit light-curves and the Rossiter-McLaughlin effect, to better constrain the orbital configuration and 
properties of planetary systems and active zones of their host stars. Tests of the code are presented to illustrate its performance 
and to validate its capability when compared with analytical models and real data. Finally, we apply SOAP-T to the active star, 
HAT-P-11, observed by the NASA Kepler space telescope and use this system to discuss the capability of this tool in analyzing 
light-curves for the cases where the transiting planet overlaps with the star's spots. 

Key words, methods: numerical- planetary system- techniques: photometry, Radial Velocity - Stellar activity 



1. Introduction 

Stellar spots and plages, combined with the rotation of 
the star, can mimic periodic Radial Velocity (hereafter 
RV) signals which ma y lead to false po s itives in the 
detec tion 
120071 : 



ay 

of a planet (IQueloz et al.l 1200 ll 



Huelamo et al. I 120081 : I Boisse et al 



Bon fils et al.l 
20091 120111) . 



Furthermore, the overlap of the planet and the active 
zones on the surface of a star can create anomalies inside 
the transit light curve that cause difficulties in constrain- 
ing the depth of transit a nd as a result, the radius of the 
planet ( Rabus et al. 2009f ). These anomalies can also lead 
to false positives in the detection of non-tr ansiting plan- 



ets by the transit timing variati on method (jOshagh et al 
20121 : ISanchis-Oieda et alll201ll ). 



The stellar spots may provide important information 
about the orbital configuration of the star's planetary sys- 
tem. One quantity that is particularly interesting is the 
angle of stellar spin axis with respect to the orbital plane 



The tool's public interface is available 



at 



http: / / www.astro.up.pt / resources / soap-t / 



of the planet. There are two methods to measure this spin- 
orbit angle: the Rossiter-McLaughlin (hereafter RM) ef- 



ing the planet's transit 


fe.e.. Ohta et alJl2005HWinn et al. 


2005; 


Narita et al.ll2007 


:lHebrard et al.l 20081: Triaud et al. 


2009, 


2010; Winn et al. 


2010: Hirano et al. 2011), and the 



study of the anomalies inside the transit light-curve due to 
an overlap of the transiting planet with an active zone of 
its host star. If long-term and continuous observations of a 
star in photometry are accessible (e.g., target stars of the 
Kepler telescope), we will be able to detect the evolution 
of these anomalies inside transits and use them to deter- 
mine the spin-orbit angle, inclination of the stellar spin 
axis, and the configuration of stellar spot s on the surface 



of th e star (longitude, lati t ude, and size) (INutzman et al 



2011: Desert et al. 


2011; 


Sanchis-Oieda et al. 


2012) 



Sanchis-Oieda fc Winnl 2011 



stellar spin axis and orbital plane, can give us a better un- 
derstanding of the configuration of the planetary system 
as well as its forma tion, dynamical evolution, and possi 



as well as its forma tion, d ynamica l evolution, and possi- 
ble migration (e.g., IQueloz et al.ll200o[ lOhta et all 12005 
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Winn et al.l 120051: iFabrvckv & Winnl 120091; iTriaud et al. 



2Qldt iMorton fc JohnsoiJl201lt Iffirano et al-ll2Q12h . 

In this paper, we present a new tool for studying sys- 
tems consisting of a rotating star with active zones and a 
transiting planet. The code, named "SOAP-T" generates 
synthetic radial velocity variations and produces light- 
curves for the system as functions of the stellar rotation 
phase. This tool is fast and as such can be utilized to simu- 
late different possibilities for different initial configuration 
in order to find the best solution for inverse problems. 

We describe SOAP-T in section 2. We will explain its 
software platform, language, input parameters, methodol- 
ogy, implementation, and output quantities. In section 3 
we present several tests of this software to demonstrate its 
capabilities. Section 4 has to do with applying SOAP-T to 
the particular cas e of HAT-P-11 and comparing its result 
with the result of Sa nchis-Oieda fc Winnl ([201 it ). We use 
this comparison as a demonstration of the capability of 
SOAP-T in analyzing the case where the planet and spot 
are overlapping. In the last section, we discuss other pos- 
sible applications of SOAP-T and conclude this study by 
presenting its future improvements. 

2. Software description 

SOAP-T has been developed based on t h e al- 
ready published code SOAP (jBoisse et al.l |2012j) 
I http://www.astro.up.pt/soap/] that simulates spots 
and plage s on the surf a ce of a rotating star. We refer the 
reader to iBoisse et al.l ( 2012f ) for more details. The code 
in SOAP-T is written in Python with calls to functions 
written in C in order to make the code faster. 



2.1. Inputs 

The code requires initial conditions which can be assigned 
by the user in an input configuration file "driver. cfg" . 

The details of the initial p arameters for the r otating 
spotted star can be found in IBoisse et al.l (j2012[ ). Here 



we explain the extra parameters that are required by 
SOAP-T. 

Stellar parameters: One new feature of SOAP-T com- 
pared to SOAP is the implementatio n of the quadratic 
limb darkening (Mandel fc Ago j |2002l ) (SOAP considered 
the linear limb darkening for the star) 



I(r) = 1 - 7 i(l -/i) -72(1 -/i) 2 



7i+72<l- (1) 



In this equation, r is the normalized radial coordinate 
on the disk of the star, I(r) is the specific intensity 
defined to reach its maximum 1 at the center of the 
star (r = 0) and its minimum at the star's limb, and 
ji = cosO = (1 — r 2 ) 1 / 2 where is the angle between the 
normal to the surface and the observer. As shown by 
equation (1), the input parameters require two coefficients 
of quadratic limb darkening (71 and 72). 



Spot parameters: The number of spots is another 
important parameter. Unlike in SOAP where the number 
of spots could be up to 4, it is possible to choose up to 10 
spots in SOAP-T. 

Planet parameters: 

Users are able to set the initial orbital configuration 
of the planet by selecting the following quantities: initial 
value of the planet's orbital period (Porb), the time of its 
periastron passage (To), its initial orbital phase (PSo), ec_ 
centricity (e), semimajor axis (a), argument of periastron 
(lj), inclination with respect to the line of sight (i), as well 
as the planet's radius (R p i an et) and the angle between the 
stellar spin axis and the normal to the orbit of the planet 
projected on the y — z plane. We denote this angle by A. 

2.2. Detail description of computation 

We restate here from lBoisse et al.l ( 2012f ) that the rotating 
spotted star is centered on a grid of grid x grid cells on 
a y-z plane. The cells take their y and z values between 
-1 and 1, normalized to the stellar radius. Each grid cell 
(PyiPz) contains a flux value and a CCF (cross correla- 
tion function). Our goal in this section is to describe 1) 
how we add a planet to the rotating spotted star, 2) how 
we compute its impact in photometry and in RV, and 3) 
how we incorporate the overlap between the stellar inho- 
mogeneities and the planet. 

To determine the effect of a planet on the RV and 
flux of a non-spotted star in each phase step, the code 
first checks whether the planet is in the foreground or 
the background with respect to the star. For that reason, 
in the first step, SOAP-T calculates the position of the 
planet for each time step. Since the output of the code is 
in the stellar rotational phase, the code uses the following 
equation to convert stellar rotation phase into time 



T=(PS + PS ) x P 7 



rot ' 



(2) 



The quantity PS is the stellar rotation phase which can be 
between and 1, PSo is the planet's orbital initial phase 
which defines the initial position of the planet, and P rot is 
the stellar rotation period in days. This equation is only 
valid for stellar rotation periods larger than the orbital 
period of the planet, which should be true for the most 
cases of transiting extrasolar planetary systems. 

In order to determine the planet's position, the code 
has to calculate A. As explained above, A is the angle be- 
tween the stellar spin axis (denoted by h*) projected on 
the y — z plane and the unit vector normal to the plane 
of the planet's orbit (shown by n p in figure 1) also pro- 
jected on the same plane. The coordinates of vector n* are 
determined using 



= Rx(^*)Ry(i*) ( 



cosi* 
sini* sin 
— sin cos?/;* 



(3) 
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where (-0*) is the longitude of stellar spin, (i*) is the 
star's inclination, and R x and R y are the rotation ma- 
trices around the x and y axes. The coordinates of the 
normal vector h p are also determined in the same fashion 
using the planet's longitude of ascending node (fi) and 
orbital inclination (i) as 



1 \ / cosi 
>, n = R x (Q)R y (i) = sin z sin fi 
/ \ — sin i cos fi 



(4) 

The angle A is then calculated using A = fi — ip*. The 
code, subsequently, aligns the projected stellar spin axis 
with the z direction by setting xp* = 180° and uses the 
longitude of the ascending node of the planet fi = 180 + 
A to calculate the coordinates of its position vector (see 
section 4 and equations 53-55 in Murray & Correia 2011 
for more details). At each stellar phase step, the position 
vector of the planet (D) is given by 



D 



sin i s'm(u + /) 
■ cos A cos(cj + /) + sin A cos i sin(uj + /) 
- sin A cos(cj + /) — cos A cos i sm(uj + /) 



In equation (5) 
_ a(l-e 2 ) 



1 + e cos / ' 



(5) 



(6) 



where r p is the distance between the stellar center and the 
center of planet, and 



f = 2 tan" 



1 + e 
1 - e 



1/2 



tan 



E 



(7) 



is the true anomaly of the planet. To calculate /, for small 
eccentricity, the code uses Kepler's equation 



E = M 



esinM - M 
1 — e cos M ' 



(8) 



where M = 2tt(T — Tp)/P or ^ is the mean-anomaly of the 
planet 's orbit (e.g. lOhta et all ( 2005 ): iMurrav fc Correia] 
( 20111 )) (see Fi gure 1, a schematic view of a planetary 
system). 

If the planet is in the foreground (D x > 0), the code 
checks if the planet is inside the stellar disk or outside by 
comparing the projected distance of the planet's center 
to the stellar center with the radius of star (|D| < R*). 
If the planet is inside, the code resolves the area of the 
grid where the planet is located and establishes an area 
around the planet where it will focus the subsequent cal- 
culations. This area is then scanned to determine whether 
each grid cell is located inside or outside the stellar disk 
(especially in the case where the planet is on the stellar 
limb), and if this grid cell belongs to the planet's disk, 



the code determines the location by calculating the dis- 
tance of the grid cell to the center of planet and com- 
paring it with planet's radius. Assuming that a grid cell 
is located inside the stellar disk and also belongs to the 
planet's disk, its CCF is modeled by a Gaussian with a 
width and amplitude given by the input parameters and 
doppler-shifted according to the projected stellar rotation 
velocity weighted by the quadratic limb-darkening law. 
This CCF value will then be removed from the total CCF 
of the rotating non-spotted star. For the same grid cell, 
the values of the flux can be calculated using only on the 
quadratic limb- darkening law. The code also removes the 
flux value of that cell from the total flux of the rotating 
spotted star. 

In the next step, the code adds a planet to the rotating 
spotted star in order to calculate the effect of both the 
planet and the spots on the RV and flux at the same time. 
In the case that the spot is not partly covered by the 
planet, the procedure is as explained below. The code 

— determines the position of the spot, 

— determines the position of the planet, 

— applies inverse rotation to the spot to bring it along 
the line of sight, 

— scans the region of the spot's position to identify the 
grid cells that belong to the spot's disk area, 

— sums up the CCF and flux of all cells inside the spot 
and removes them from the total flux and CCF of the 
non spotted-star [see figure 4 of Boisse et al. (2012) for 
the inverse rotation and also for more details on spots 
calculation] , 

— scans the region of the planet's position to determine 
the grid cells that are inside the disk of the planet, 

— calculates the sum of the CCF and flux for those grid 
cells inside the transiting planet and removes them 
from the total flux and CCF of the spotted star. 

The complexity arises from the overlapping of the tran- 
siting planet with a spot on the surface of the star. In 
this case, the code considers those parts of the spot whose 
distances to the center of the planet are smaller than the 
planet's radius as the overlapped area. Those points will 
not be scanned during the spot scanning process. As a 
consequence, they will skip the rest of procedure for CCF 
and flux of spot and removal from total CCF and flux. 
This procedure is able to produce those positive "bump" 
anomalies inside the transit light curve that are due to the 
reduction in the loss of light. The schematic view of the 
system which is simulated by SOAP-T is shown in Figure 
2. 

2.3. Outputs 

The code returns the results of its simulation into the 
output file "output.dat". The output file contains four 
columns that represent stellar phase, flux, RV and BIS 
(the bisector span), respectively. These quantities can be 
plotted directly by SOAP-T as functions of the stellar ro- 
tation phase. Similar to SOAP, the stellar phase step can 



4 



M. Oshagh et al.: SOAP-T. 



r t 




Ascending node 



Fig. 1. Left: Schematic configuration of the stellar spin axis and planetary orbital plane. The vector h p is the normal 
unit vector of the planetary orbital plane, the vector n* is the stellar spin axis, and i* is the stellar inclination. We 
denote the angle between the projected stellar spin axis and the projected normal unit vector of the planetary orbital 
plane on the y — z plane with A (not shown here). Right: Schematic view of planetary orbit from above the orbital 
plane, the star is the gray circle and the planet is the black full circle. See Section 2 for more details on the symbols. 



be chosen in two ways, first by selecting the constant frac- 
tion of the phase between to 1, or by uploading the 
wanted phase values in the phJn file. The flux values cor- 
respond to the photometry of the system (the total flux of 
system) , and they are normalized to the maximum value of 
total flux during the stellar rotation phase. The RV values 
are obtained by fitting a Gaussian func t ion to each total 
CCF [see more details in iBoisse et al.l (j2Q12f )]. We note 
here that each total CCF is also normalized to its maxi- 



mal value. We refer the reader to lQueloz et al.l (|200l[ ) for 
details on BIS calculations. 



3. Tests 

To check the capability of SOAP-T and the validity of its 
result, we preformed several tests to compare its outcome 
with theoretical models and real observations. Tests re- 
lated to the effects of spots on the RV and flux s ignals 
have been presented and explained in lBoisse et al. ( 2011 



2Q12h . In this section, we test the validity of the results for 
the photometry and RV with a transiting planet around 
both a non-spotted and a spotted star. 

To begin with, we consider a system consisting of a 
non-spotted star with a transiting planet. The purpose 
of this test is to determine whether the code can repro- 
duce the light curve of a star with a t ransiting planet 



as obtained from the theoretical model of Mandel fc Agol 



(|2QQ2f ) . We select various initial conditions including a star 
without any limb darkening, or with linear or quadratic 
limb darkening. As shown in Figure 3, the results obtained 
from SOAP-T are in strong agreement with the results 
obtained from the theoretical model base d on the formal- 
ism presented by iMandel fc A^ol (|2002l) . The latter im- 
plies that SOAP-T can be used reliably to characterize 
the light-curve of a transiting planetary system. 



In a second test, we examined the RVs obtained us- 
ing SOAP-T during the transit of a planet in front of 
a non-spotted star. We applied SOAP-T to the system 
of WASP-3 which is known to harbor a planet with a 
misalignment between its orbital plane and stellar spin 
axis (see the parameters of WASP-3 system in the Table 
1). The RV observatio n data of WASP-3 du ring its tran- 
sits were taken from ISimpson et al.l ( 2010f ) (it was ob- 
served with the SOPHIE spectrograph at Haute-Provence 
Observatory). We fitted the RVs produced by SOAP-T 
to the o bservational data and compared the results with 
those of lSimpson et al.l (|2010l ) . We only allowed stellar ro- 
tation velocity (v sini) and misalignment angle (A) to vary 
as free parameters, and kept all other parameters constant 
and equal to their values reported bv lGibson et al.l ( 2008 ) 
and Simpson et al. ( 2010f ) (see Table 1). 

The best fit result obtained using SOAP-T corresponds 
to A = 20.0° ± 3.3° and v sin i = 13.31 ± 0.45 kms" 1 with 
Xred = 0-9308. We obtained the error bars of these mea- 



surements using the bootstrap method ( Wall fc Jenkinsl 



2003). The value of the misalign ment angle (A) i s con - 



sistent with the value reported by ISimpson et al.l ( 2010f ) 



(13°±7o) within la. Note that the error bars obtained 
with SO AP-T are more than twice smaller than those de- 
rived by ISimpson et al. (2010). We assume that this could 
be due to smaller number of free parameters in our fit- 
ting procedure, otherwise it might be also due to the fact 
that SOAP-T has been developed to reproduce exactly the 
analysis routine of RV measurements and is thus closer to 
the observations when compared to the resu lt of the ana- 
lytica l formula used by ISimpson et al. (l2010h dBoue et al 



20121 ). Indeed, we obtained also a better agreement be- 



tween the projected stellar rotation velocity (vsini) and 
the value obtained by spectroscopic broad ening measure- 
ments (13.4 ± l.bkms~ l ). For comparison. ISimpson et aT 
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JD- mid transit 



Fig. 2. Schematic view of stellar spot and the transiting 
planet in simulations by SOAP-T. 



(2010) obtained vsini = 19.6 ±2!? kms 1 . Figure 4 shows 
the result of the best SOAP-T RV's fit and the observation 
of RM effect of WASP-3. 

Finally, since there is no theoretical model for the case 
where a transiting planet and star spots overlap, to eval- 
uate the performance of SOAP-T in this case, we com- 
pared the photometric result of SOAP-T with the results 
reported in the literature for the real observation of HAT- 
P-ll by the Kepler space telescope. We will explain this 
test and its applications in the next section. 

4. Application of SOAP-T to HAT-P-11 

HAT-P-11 is a 6.5 Gyr bright (V=9.6) K4 star with a mass 
of 0.81 Mq and radius of 0.75 R@. At a distance of 38 
pc, this star hosts a transiting Neptune-sized planet with a 
mass of 0.081 Mj up and radius of 0.422 Rj up (IBakos et al 



20101 ). HAT-P-llb has a period of 4.887 days and revolves 
around its central star in an orbit with a semimajor axis 
of 0.053 AU and an eccentricity of e = 0.198. More details 
of the parameters of this planet and its host star can be 
found in Table 2. 

Studies of the HAT-P-11 planetary system have shown 
extreme misalignment of the projected angle between the 
stellar spin axis and the pl anet's orbital pla ne (A = 
103°±?go) through RM effect ([Winn et al.ll2010l ). 

HAT-P-11 (KOI-3) was used as one of the calibration 
target stars for the Kepler space telescope and has been 
observed throughout Kepler's mission in both short and 
long cadences. The high precision photometry of this star 
has revealed some features inside its transit light curve 
and has also shown a large modulation outside of transits. 
The anomalies inside the transit light curve of HAT-P- 
11 can be attributed to the overlapping of the transit- 
i ng planet with active zones on the surface of the star 
(jSanchis-Oieda fc Winnll201ll ). The recurrence (or not) of 



Fig. 3. Comparing the photometric results of SOAP-T 
and the theor etical model of a transi ting planet over a non- 
spotted star ( Mandel fc Agoll 12002 ). The cyan, red, and 
yellow lines show SOAP-T's results for a star without limb 
darkening, a star with linear limb darkening (71 = 0.6), 
and a star with quadratic limb darkening (7I = 0.29 and 
72 = 0.34), respectively. The dash-dotted line, dashed 
line, and the dotted lin e show the correspond ing results 
using the mechanism bv lMandelfc Agoll (120021) . 




0.4B 0.50 0.52 

Orbital phase 



Fig. 4. RV observations of WASP-3 during transit of 
WASP-3b (RM effect) minus K epler ian motion, which are 
blue dots ( Simpson et al.llioioL overplotted with the best 
fit of RV signals obtained from SOAP-T, which is red line. 



these anomalies can be us ed to put an upper limit on the 
stellar obliquity. Recently ISanchis-Oieda fc Winnl ([201 lh 
carried out a purely photometric study of this system 
based on the observation of photometric anomalies inside 
the transit light curve using the Kepler public data Q0- 
Q2 and showed that no recurrence of anomalies exists in 
two closely spaced transits in the Kepler data of HAT-P- 
11. They considered this result as evidence of misalign- 
ment of the stellar spin axis with respect to orbital plane. 
Fortunately, these authors were able to detect a feature 
with two peaks in the folded light curve of 26 transits of 



6 M. Oshagh et al.: SOAP-T. 

Table 1. The stellar and planet's parameters of the WASP-3 system ( Gibson et al.ll2QQ8t ISimpson et al.ll2010f ). 



Parameter 


Value 


Uncertainty 


Stellar parameters 








1.31 


i 0.05 
^0.07 


Linear limb darkening coefficient 


0.69 




Stellar inclination (deg) 


90 




Spectroscopic stellar rotation velocity v sin i(kms~ 1 ) 


13.4 


±1.5 


Planet parameters 




_l_0.0014 
^0.0013 


Planet to star radius ratio (R p /R*) 


0.1013 


Period (days) 


1.846835 


±0.000002 


Eccentricity 







Scaled Semimajor axis (a/R*) 


5.173 


4_0.246 
^0.162 


Orbital inclination (deg) 


84.93 


,1.32 
^0.78 


Projected spin-orbit misalignment angle(deg) 


13 


,9 
±7 
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Fig. 5. Comparing the best fit model to the transit pho- 
tometry of HAT-P-11 using SOAP-T and the Kepler ob- 
servation of HAT-P-11 for one period of stellar rotation. 
Red (dash-dotted) line shows SOAP-T's photometric re- 
sult for an edge-on solution and green(dashed) line corre- 
spond to HAT-P-11 observation. The blue (solid) line in 
the bottom panel shows the residual. 

this star which could be interpreted as evidence for the ex- 
istence of two long lived spot belt regions on the star. By 
using a simple geometric model, these authors arrived at 
two solutions for the stellar obliquity; an edge -on solution 



that i s in good agreement with the result of IWinn et al 
(2010) through the RM study, and an alternative solu- 
tion in which the star is seen almost po le-on. The method 
used by ISanchis-Ojeda fc Winnl ( 2011 ) also allowed these 
authors to put some constraints on the position (latitude) 
and size of the active zones on the surface of star (see 
details of parameters for two solutions in Table 3). 

Here we present the results of the application of SOAP- 
T to HAT-P-11 in order to reproduce the light curve of 
this star over one stellar rotation. The rotation period of 
this star contains six transits. 

Since the inside transit anomalies carry more informa- 
tion about the configuration of a planetary system, we 
decided to give more weight to these points in our study. 
We, therefore, performed a regular sampling to reduce the 



number of points outside of the transit light curve. One 
point was taken from each 300 points in the light curve out 
side the transit. This process reduced the running time of 
the code significantly. 



4.1. HAT-P-11 Edge-on solution 



In this section, we consider the system to be edged-on and 
choose all initial parameters of the star, planet, and spots 
(except the number of spots and their long i tude) from 
the values reported bv lSanchis-Ojeda fc Winn ( 2011 ). The 
spot's brightness are fixed to zero. The details of initial 
conditions are listed in Table 3. To reproduce all the fea- 
tures in the light curve, both inside and outside the tran- 
sit, we limited the la titude and the size of the sp ots to 
the range reported bv lSanchis-Ojeda fc Winnl ( 2011 ). The 
spots's radius are equal to the half-width of the active 
zone which correspond to R spo t = 0.08i2*). However, we 
allowed the number of spots and their longitudes to vary as 
free parameters between 0-10 and 0° — 360°, respectively. 
The best fit to the observations (xled = 7.356) was ob- 
tained with eight spots on the surface of the star. Figure 5 
shows the best result of SOAP-T photometry overplotted 
with the real observation of HAT-P-11 during one stellar 
period. Figure 6 shows more details of the inside of each 
transit and the results obtained from SOAP-T after nor- 
malizing both flux values to one. As shown in both figures, 
SOAP-T could manage to reproduce all the inside tran- 
sit anomalies and also the out side transit variation. We 



note here that the results of lSanchis-Oieda fc Winnl (|201l[ ) 
were obtained by considering the inside transit anomalies. 
These authors did not consider the outside features of the 
light curve. However, using SOAP-T we were able to study 
the entire light-curve, both inside transit and outside, at 
the same time. We n ote that here we do no t inten d to 
confirm the results of Sanchis-Ojeda fc Winnl ( 2011 ). We 
only use the similarities between the results obtained from 
SOAP-T photometry (e.g. position and size of spots, stel- 
lar inclination and spin-orbit misalignment angle) and the 
results by these authors as a proof of the validity of SOAP- 
T's capability in simulating the case of spot and planet 
overlapping. 
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Tabl e 2. The stellar and planet's parameters of the HAT-P-11 system ( Bakos et al.ll201Q[ ISanchis-Oieda & Winn 
201 lh . 



Parameter 



Value 



Uncertainty 



Star parameters 

M*(M Q ) 0.81 

R*(Rq) 0.75 

Stellar rotation period (days) 30.5 

Age (Gyr) 6.5 

Distance (pc) 38.0 

Linear limb darkening coefficient 0.599 

Quadratic limb darkening coefficient 0.073 

Planet parameters 

M P (M Jup ) 0.081 

Planet to star radius ratio (R p /R*) 0.05862 

Orbital period (days) 4.8878049 

Scaled Semimajor axis (a/R*) 15.6 

Eccentricity 0.198 



i 0.02 
^0.03 

±0.02 

,4.1 
±3.2 

.5.9 
=C 4.1 

±1.3 
±0.015 
±0.016 



±0.009 
±0.00026 
±0.0000013 
±1.5 
±0.046 



Table 3. Parameters of the HAT-P-11 system for different solutions reported bv lSanchis-Ojeda fc Winn ( 2011 ) 



Parameter 


Value 


Uncertainty 


Edge-on solution 






Projected spin-orbit misalignment angle A (deg) 


106 


+ 15/-12 


Stellar inclination i s (deg) 


80 


±4/-3 


Latitude of active zone (deg) 


19.7 


+1.5/-2.2 


Half-width of active zone (deg) 


4.8 


+ 1.5/-1.8 


Pole-on solution 






Projected spin-orbit misalignment angle A (deg) 


121 


+24/-21 


Stellar inclination i s (deg) 


168 


±2/-5 


Latitude of active zone (deg) 


67 


±2/-4 


Half-width of active zone (deg) 


4.5 


+ 1.6/-1.9 




'J 





Fig. 6. Zoom on the six transits of HAT-P-11 in the figure 5. The flux values are normalized to one for both observation 
and SOAP-T. Red dot represent the SOAP-T photometric result and green triangles are observed flux of HAT-P-11. 
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4.2. HAT-P-11 pole-on solution 

In this section, w e examine the pole-on model f or the sys- 
tem proposed by ISanchis-Oieda fc Winnl (|201l[ ) . We note 
that here we explore which configuration of spot on the 
surface of the star would be able to generate the out of 
transits variation and we did not attempt to reproduce 
the inside transit features. We generated the light curve 
of HAT-P-11 using SOAP-T photometry and the parame- 
ters of the pole-on model given by these authors, (see the 
details of paramete rs in Table 3). Since the pol e-on so- 
lution presented bv lSanchis-Ojeda fc Winn ( 201 lh is only 
nearly pole-on (i.e., the angle of misalignment is ~ 168°), 
some spots located close to the equator of the star can 
produce large variations outside the transit by disappear- 
ing from the view for a small part of the rotation period. 
We investigated the effect of different possible sizes of the 
spots that were located in different latitudes. Our analy- 
sis showed that the optimized result, whose outside-transit 
profile is similar to that from the observation, has to have 
a spot on 20° or 60° latitude with a size equal to 0.2R* 
(see Figure 7). Since the brightness of spot is fixed to zero, 
it is the minimal size of the spot. For a star with a sim- 
ilar spectral type as HAT-P-11 (i.e., K), this is a very 
large spot. Although s pots with very lar ge sizes have been 
detected on K stars ( Strassmeierl [l999) , this situation is 
unlikely for HAT-P-11. 



5. Conclusion and Perspective 

In this paper, we introduced a new software package that 
can be used to produce the photometric and RV signals 
of a system consisting of a rotating spotted star and a 
transiting planet where the spots and planet overlap. We 
tested the capability of our code by comparing its results 
with theoretical models and the results of observations. 
We used the HAT-P-11 system and demonstrated that 
when using the edge-on model of ISanchis-Oieda fc Winr] 
()2Q11[ ). SOAP-T is capable of reproducing the same fea- 
tures as obtained from observations for the outside of tran- 
sits as well as the anomalies inside the transits. We also 
showed that the reconstruction of the large modulation of 
th e outside transit of HAT-P-11 with the pole-on model 
of Sanchis-Ojeda fc Winnl ( 2011 ) requires an assumption 
on the size of spot on the star which may not be realistic. 

In using SOAP-T, the number of free parameters of a 
system that can be chosen by the user is large. In gen- 
eral, for a system consisting of a rotating star with one 
spot and a planet transiting, there are 18 free parame- 
ters. As a result, for a system where parameters are not 
constrained (e.g., by observation), the finding of the best 
solution that fits real observation becomes a complicated 
task. The minimization of the reduced chi-squared (the 
indicator of goodness of the fit) in this case requires the 
exploration of a large parameter space. A simple griding 
of this space turns the calculations into a time consuming 
process which may result in the finding of local minima in- 
stead of a global minimum. Efforts are currently underway 



to investigate different optimization strategies which can 
lead to faster and more secure determination of a global 
minimum. 
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Fig. 7. Comparing SOAP-T' s photometric result with the observation of HAT-P-11 for the pole-on model of 
Sanchis-Oieda fc Winnl (|201lf ). Each panel shows different values for the latitude of one spot with the size of 0.2 
radius of star. Green dot is the SOAP-T photometric result and red circle is the observation of HAT-P-11. 
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